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A brief outline of Dark Matter detection experiments using liquid Ar- 
gon technology is presented. The Pulse Shape background discrimination 
method (PSD) is described and the example of its use in 2.3 1 R& D detec- 
tor is given. Methods of calculating sensitivity of a Dark Matter detector 
are discussed and used to estimate the possible improvement of sensitivity 
after introduction of isotopically depleted liquid Argon. 
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1. Dark Matter Detection 

1.1. Experimental evidence and candidates 

Since the beginning of the 20th century the number of experimental 
facts suggesting the insufficiency of using solely luminous matter in the 
description of the Universe has been growing. Several measurements have 
shown that the baryonic matter provides only the 4.5% of the mass - energy 
of the Universe. Of the remainder: about 73% is associated with Dark 
Energy (Cosmological Constant) and over 22% with Dark Matter 

One of the prominent observations suggesting the existence of DM was 
the measurement of the shape of the rotation curves of galaxies. The 
stars rotate around the center of the galaxy with such velocities as if the 
mass was distributed in an uniform way, and not in the central bulge and 
the surrounding disk. The observation of so called "polar ring galaxies" [4J 
(with additional stars rotating perpendicularly to the main disk) leads to a 
similar conclusion. 
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Fig. 1. One of the evidences for the Dark Matter existence: the Buhet Cluster 




Fig. 2. Another important evidence for the Dark Matter existence: the map of the 
CMB distribution from the WMAP experiment [3]. 

Similar conclusions can be drawn from the observations of objects such 
as the Bullet Cluster (Figure [l] right) where the movement of the luminous 
parts does not trace that of the main mass component of the colliding clus- 
ters as measured by weak lensing. The measurements of Cosmic Microwave 
Background Radiation (CMB), Figure [2j also support the hypothesis of 
Dark Matter. The anisotropics in the CMB can help calculate of the per- 
centage of the baryonic matter in the Universe and the share of the Dark 
Matter in its total mass. 

Today, most astrophysicists agree that Dark Matter exists - it is its na- 
ture that is the question. The number of candidate particles is great. These 
can be divided into warm DM (the relativistic particles which freezed out at 
the early stages of the Universe formation), which will not be discussed here, 
and the cold DM which includes particles with non - relativistic velocities 
and yet not known form. The last is included in the preferred cosmolog- 
ical model - ACDM. The most popular and promising DM candidate is a 
WIMP (a Weakly Interacting Massive Particle). Its allowed mass range ex- 
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tends from tens of GeV to above 1 TeV [6J . The excluded cross sections are 
of over 10"^"^ cm^ for spin independent interactions and over 10~^^ cm^ for 
spin - dependent interactions [5j. As the velocities are non relativistic, the 
expected energies of recoiling nuclei in a detector are of the order of tens 
of keV [7J. Candidates for WIMPs may be supersymmetric particles (e.g 
neutralino). The expected event rate is very low - 0.02 events/day/kg of 
detector mass. The issue of Dark Matter detection has to be dealt with by 
employing efficient detecting techniques and and reducing the background 
to unprecedented levels. 

1.2. Overview of current detection techniques 

The methods of Dark Matter detection can be divided into indirect 
and direct. The former uses products of annihilation or decay of the DM 
particles and is not discussed here, whereas the latter is connected with 
the observation of DM interactions in a detector. The transfered energy 
can be detected as ionization of the detector medium, heat in the crystal 
lattice of the detector or scintillation (examples of experiments are given in 
Figure [s]). Typically, a combination of two of the mentioned signals is used 
to obtain a better background discrimination. Cryogenic detectors using 
Ge or Si crystals (CDMS) detect ionization and heat, detectors based on 
CaW04 use the scintillation signal and heat and liquid noble gas detectors 
- ionization (when working in double phase mode) and scintillation. 




CaWO,: CRESST, 
ROSEBUD 



Fig. 3. Classification of experiments according to their detection techniques [8] 



2. Liquid Argon as a Dark Matter detector 

The analysis presented in this work uses some of the methods developed 
in [9] and expands on the results presented in [10] and as such, will focus on 
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liquid Argon detectors. A typical double-phase Argon DM detector consists 
of a dewar filled with liquid argon with photomultipliers observing a fiducial 
volume. Over the liquid, there is a layer of gaseous Ar. 




Fig. 4. A scheme of a LAr DM detector (based on [12]). 



The detector can work in single phase mode, registering only the light 
signal produced in the liquid (SI). In this detection mode the only sig- 
nal which is measured is the scintillation. When a DM particle hits Argon 
nucleus, it may become excited or ionized and form a dimer Ar2. The relax- 
ation of the excited dimer may be of two forms - singlet (fast signal) with 
decay time about 7 ns and triplet (slow signal) with decay time 1500 ns. 
Such a significant difference (absent in Xenon, also used in DM detection) of- 
fers a background discrimination possibility because electrons produce slow 
and fast components in a different proportion than neutrons (and WIMPs). 
This results in different signal shapes. When the double phase mode is 
used, the ionization caused by the recoiling nuclei is extracted to the gas 
layer above the liquid by the electrodes situated below and over the surface 
to cause the secondary signal (S2) proportional to the ionization (Figure [4] 

Using SI and S2, the difference in the ratio of relative amplitudes of 
signals resulting from electron-like events (main background) and neutron 
(or WIMP)-like events is an additional discrimination method, resulting in 
a better rejection of background events (Figure [5]). 



2.1. Background suppression 

The background in DM detectors can be divided into two subclasses: 
external and intrinsic. Both are caused mainly by 7 rays and neutrons. 
The effects of the background resulting from the cosmic rays (dominating 
during tests performed on the earth surface) and the natural radioactivity 
of the rocks (in underground laboratories), can be reduced by precise mea- 
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Fig. 5. A difference of the SI and S2 signals for an electron - like event (left) and 
the neutron - like event (right) [9]. 



surements and simulations leading to determination of their signature as 
well as by using shielding. 

The intrinsic background connected with the construction of the detector 
can be diminished by using radiopurified materials, but the most important 
part of it is related to the radioactive isotopes ^^Ar and ^^Ar present in at- 
mospheric Argon along with non radioactive "^^Ar. Produced in interaction 
of cosmic rays or due to the neutron capture, these isotopes contribute to 
the background due to /3~ decays [13j. The activity of ^^Ar dominates and 
was measured to be 1.01 it 0.02(stat) ± 0.08(syst) Bq for 1 kg of liquid 
Argon [14j. 

To reduce background effects, in double phase detection mode the S1/S2 
discrimination method can be applied. In single phase measurements, like 
presented in this work, the Pulse Shape Discrimination has to be applied. 
We have used the FPrompt method described below, as it offers good dis- 
crimination power and its implementation is not difficult. The studies which 
led to the choice of FPrompt are described in [lOj. The parameter used to 
separate background and signal is calculated according to the formula: 

Fp = ^ = % ^ (2.1) 

where Ti and Tp are the limits of the time window for signal V(i) and Tp^ 
is the time which ensures the best separation of two populations - 100 ns 
for analysis of 2011 data and 120 ns for the data from 2009 [lOj. 



2.2. Detector upgrade performed in 2010/2011 

The WArP 2.3 liter prototype detector (Figure [7| left) was the first 
to report DM search results in liquid Argon [15j and has been since used 
for R &D purposes. In 2009 the WArP collaboration measured the n-like 
ev./e~-like ev. separation by irradiating a chamber with an Am/Be source. 
In the FPrompt distributions calculated using this data, an intermediate 
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population was found. It was connected with the inelastic scattering of 
neutrons and included in the model which was then fitted to the spectrum. 
Indeed, such scattering dominates for energies of the particles emitted by 
the Am/Be source (Figure [g]). The data suffered from a low light yield and 
due to some DAQ problems it was not possible to subtract the background. 
After a detector upgrade, we decided to repeat this measurement profiting 
from a four times higher LY. This upgrade, performed in 2010, is described in 
detail in [16j and [17j. New photomultipliers with higher quantum efficiency 
were installed, which led to a higher light yield (6.1 phel/keV compared to 
1.52 phel/keV in 2009 runs). The consequence of this was an insight into 
a region of lower recoil energies. The new DAQ electronics with a broader 
dynamic range, causing less saturated events, made the exploration of higher 
energies possible (Figure [t] right). 









1 












> 


















c 






V o.z 

0.15 
0.05 






c 


2 4 t .5 10 




Energy (MeV) 




Fig. 6. The spectrum of the Am/Be source (left) [TP and references therein and 
the cross section for neutron interaction on argon nuclei (right) [9^ and references 
therein. 




Fig. 7. The the detector scheme (left) p!8l and a difference of the energy ranges in 
2009 and 2011 data (right). 
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2.3. Analysis of 2011 data and reanalysis of 2009 data 

In 2011, runs were performed to test the new detector setup. The Am/Be 
neutron source was used to imitate the WIMP interactions and the ^^^Am 
to perform detector cahbration and to estimate the hght yield. The results 
from 2009 were used as a reference. The FPrompt spectrum was fitted with 
a sum of two Gaussian functions (for gamma and neutron population) and 
a convolution of an exponential and Gaussian function (for the intermediate 
- inelastic neutron interaction population): 



Gn^iGilS)Ei)^G^ (2.2) 

where the Gn and are the Gaussian distributions for a neutron and 
gamma population respectively and Gi Ei is the convolution. One of the 
purposes of this analysis was to test whether the intermediate population 
would be observed (one of the energy bins is presented in Figure [s]). 
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Fig. 8. The FPrompt spectrum with an intermediate population visible ( this anal- 
ysis results). 
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The positions of the gamma and neutron populations resulting from 
these fits are compared with the previous measurements and published re- 
sults from other groups ([lOj, [19j) in Figure [9] 

In order to understand the small discrepancy between the new data and 
the 2009 results, a re-analysis of the 2009 data with the use of new scripts 
and reconstruction procedures was performed. The preliminary results are 



presented in Figures [TO] and [TTJ Figure [TT] presents a comparison of the 
number of elastic neutron interactions per unit of energy to the number 
expected from Monte Carlo simulations. 



8 



v43p1509 printed on October 10, 2012 





2011 analysis gammas i — i — i 




2011 analysis neutrons i — ^. — i 




Ijppincott et al. gammas > — — > 




Uppincott et aJ. neutrons ■ — b — i 




2009 analysis gammas ■ ■ ■ - 




2009 analysis neutrons ' — * — 











200 400 600 800 1000 1200 1400 



energy [keV] 

Fig. 9. Comparison of the FPrompt spectrum for the 2011 and 2009 WArP data 
[9] and the results published by Lippincott et. al |l9j . 



reanalyzed neutrons 




Fig. 10. The results of the re-analysis of the 2009 data with the new reconstruction. 



3. Predicted Sensitivity for a Depleted Argon Detector 

3.1. Calculating the Sensitivity of a Dark Matter Detector 

Even if a DM detector sees no signal, it is still possible to set the ex- 
clusion limits on the allowed a vs WIMP mass parameter space. Typically, 
the model accounting for the Earth galaxy motion, the velocity distribution 
of WIMPs in the galactic halo, seasonal DM flux change and detector effi- 
ciency is applied, as described e.g in [7J. The resulting function predicting 
the number of observed events as a function of recoil energy is integrated 
taking the detector threshold into account. Repeating this procedure for all 
WIMP masses and taking into account the detector exposure leads to the 
sensitivity plot. The tools using this approach were prepared and success- 
fully tested on published data. Another way to calculate detector sensitivity 
was proposed by S. Yellin [20] and is based on the distribution of events as 



v43pl509 printed on October 10, 2012 



9 



Real 6t Simulated Neutron Population | 



1 




I 




• F^eal NnjtrQn Papulation 










0.8 
0.6 




"l. 






0.4 











0.2 







■ 1 








< < . < 


■ . • 

' " : f 


5 


1C 


150 200 250 



Energy (phe) 



Fig. 11. Comparison of the neutron population share in the fitted FPrompt spec- 
trum. The black dots is the 2009 data, the red - the MC simulation (figure from 
[To]) and the green - the 2011 data without background subtraction applied added 
to the plot. 



a function of their energy. This was also tested but not used in presented 
analysis. 



3.2. Depleted Argon 

One of the crucial parameters in calculating the sensitivity of a DM 
detector is its exposure. Although the sensitivity should increase with the 
exposure, it is limited due to the presence of radioactive ^^Ar isotope. Its 
abundance causes the pile - up of the signal and leakage of the gamma 
population (background) into the signal region. The isotopically depleted 
Argon seems to be a promising solution. The abundance of the ^^Ar can be 
reduced by depletion in centrifuges but costs of such production method are 
high. The other approach is to search for underground Argon, which may 
be found e.g in reservoirs of Helium and then distilled. Using this method, 
a depletion factor of 25 has been confirmed [2T] . 



3.3. Predicted sensitivity for a Depleted Argon detector 

The method used in this work is based on the known ^^Ar energy spec- 
trum and radioactivity. The maximum exposure before the background 
events start to "leak" into the conservatively chosen signal region is calcu- 
lated. To do this, the events are randomly chosen from the analysis energy 
range and according to the known parameters of the Gaussian function in 
the corresponding bin the number of entries in the signal window is esti- 
mated. This, thanks to the known radioactivity of ^^Ar, helps to estimate 
the possible exposure (Figure 12). 
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Fig. 12. The method used to estimate exposure (left) and the ^^Ar spectrum with 
the analysis range marked in green (right). 



The calculations of the predicted exposure have been performed assum- 
ing ^^Ar to be the only background. The results from XENON (2011) were 
used as the reference and are presented in Figure [l3] ([22], 122 )• The case 
with the S2/S1 discrimination applied (with conservative assumption of the 
background rejection power of 10^) is displayed with a green curve and the 
case with a depleted Argon (by factor 25) - with a blue curve. 




10^ 10^ 10^ 

WIMP Mass [GeV/c^] 



Fig. 13. The comparison of sensitivity obtainable in a liquid argon detector with 
(blue) and without (green) the depleted Argon. Results from XENON 100(2011) 
[23 were used as the reference (red curve) and plotted using tool [22 . Other curves 
are obtained using tools prepared for this analysis. 
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4. Conclusions 



The preliminary results for the analysis of the data obtained with the 
upgraded detector setup were presented. The insight into a wider energy 
range has been reported. The comparison with the data from 2009 has been 
performed. The influence of the possible use of the isotopically depleted 
Argon has been demonstrated and the significant possible improvement of 
sensitivity has been seen. 
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